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Abstract. The study was conducted in the Polish (Roztoczanski National Park) and Ukrainian (Rava-Rus’ka Landscape 
Reserve and Yavorivskyi National Park) parts of the Roztocze region. In each of these locations three research areas 
were established in Scots pine (Pinus sylvestris L.) stands under similar ecological conditions. The purpose of this 
study was to carry out a survey of possible scenarios for pine stand dynamics in the Polish and Ukrainian parts of 
Roztocze using the FORKOME model. A control scenario was compared with four other climate change scenarios 
(warm dry and warm humid; cold dry and cold humid) for a period of covering the next 100 years. Using the control 
scenario, the FORKOME model predicted that for the next 100 years pine stands will dominate in terms of biomass 
and number of trees. The warm-dry and warm-humid climate scenarios resulted in slightly reduced biomass of pine 
stands. However pine would still maintain its dominance, although with a noticeable increase in beech and fir biomass. 
Nevertheless, in term of the number of trees during the second half of the simulation, it is beech and fir that dominate 
stand 1 in Roztoczanski National Park. Under the climate cooling scenario (cold dry and cold humid), the biomass of 
pine and spruce would increase during the next 100 years. Pine trees that would dominate in terms of their numbers, 
although the number of spruce individuals also tends to increase. The results presented in the paper indicate that the 
FORKOME model is very useful when investigating different climate changes scenarios in the Roztocze region. 
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1. Introduction 

Dynamics of changes in pine stands and defining 
their directions regarding possible climate change 
conditions in Roztocze are important issues in forest 
management, both from the theoretical and practical 
point of view. 

There is lack of data in literature regarding computer 
simulation of the dynamics of pine stands in Roztocze. 
Studies focusing on general analyses of pine stands 
mention that the area of Leucobiyo-Pinetum decreases 
as a result of intensive increase in the number of beech 
and fir. Moreover, many substitute forest communities 
were formed as a result of previous forest management 
geared towards pine production (Izdebski et al. 1992). 



In addition, progressing eutrophication of forest stands 
causes decreasing of Leucobryo-Pinetum communities 
and increasing range of mixed forest (Izdebski et al. 
1992). 

Other studies (Maciejewski 2009) confirm 
increasing deciduous species at the cost of decreasing 
coniferous species in the Roztoczanski Narodowy Park 
(Roztoczanski National Park - RNP). 

Catastrophic weather phenomena, for example frost 
and drought, wind-fallen trees, snow-fallen trees, and 
anthropogenic factors, are some of the reasons of the 
aforementioned changes. 

Significant influence of climate fluctuation and 
industrial pollution may be noticed. It directly and 
indirectly influences eutrophication of forest habitats, 
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and enables expansion of termophilic species (beech, 
fir) that intensively are shading the soil. It is important 
to test future dynamics of pine stands regarding climate 
changes, especially to test prognosis of pine and other 
species which share the forest stands by applying 
computer modelling. 

The FORKOME model (Kozak et al. 2012) was 
applied to predict changes of pine stands. The model 
is similar to the JABOWA model (Botkin et al. 1972), 
FORET model (Shugart 1984) and the FORLAS model 
(Brzeziecki, Zajqczkowski 2006). 

The aim of the study was to carry out prognosis of 
possible changes of pine stands in Polish and Ukrainian 
parts of Roztocze during the next 100 years. Different 
scenarios of changes were included: a control scenario, 
warm-dry climate and warm-humid climate, cold- 
dry climate and cold-humid climate. We carried out 
a simulation of dynamics of pine stands with the use 
of FORKOME model on the basis of 3 research areas 
of pine stands grown on podzol soil. Influence of 
climate change conditions (changes of total effective 
temperatures and yearly total precipitation), decreasing 
process of pine and rebuilding of other species in this 
place (beech, fir, spruce, and hornbeam). Climate 
scenarios were to verify a thesis of pine maintained 
on researched areas in Roztocze under climate change 
variants. Aggregation of climate factors approach (Kahn 
1994) was included in the analyses. 

2. Materials and methods 

The described research was performed between 
2009 and 2011 in the Roztoczanski National Park, 
Poland (RNP), in the Poland and also in the Yavorivskyi 
National Park, Ukraine. 

Three research areas were established at each of 
the aforementioned locations in forests with similar 
ecological conditions. 

One research area from each location was selected 
for presentation herein. Proper podzol soils, poor in 
minerals prevail at all the locations with the Leucobryo- 
Pinetum plant community, which is the most common 
type of forest association in the whole Roztocze 
region. Forest stands characterised by similar species 
composition (10 pine + beech, fir, spruce) and age (one 
birch which vanished in the first years of prognosis was 
observed in the first area) prevailed in all research areas. 
The research areas were located in different parts of the 
Roztocze region. The first plot - ‘Pine 1’ is located in 



the forest district no. 189f in the Roztoczanski National 
Park (50°36’N, 22°57’E), Poland, the second plot ‘Pine 
2' is in the forest district no. 66/3, in the Rava-Rus’ka 
Landscape Reserve (50°13’N, 23°37’E), Ukraine, while 
the third plot ‘Pine 3’ - is in the forest district no. 28/12, 
Mlynkivskyi forest range of the Yavorivskyi National 
Park, Ukraine (50°03’N, 23°49’E). 

Scots pine ( Pinus sylvestris L.) dominated and 
common beech ( Fagus sylvatica L.), fir ( Abies alba 
Mill.) as well as Norway spruce ( Picea abies L. Karst.) 
were present in all the areas. Canopy cover in layer 
A did not exceed 70-80%. Due to the death of one or 
a few larger Scots pines, the cover was not distributed 
evenly. Groundcover consisted predominantly of 
common bilberry ( Vaccinium myrtillus L.), grasses and, 
locally, mosses. 

Pine, the prevailing tree species, grew in dry habitats 
and formed sunny and bright stands. Spruce, beech and 
fir were present in the bush layer. For pine, the mean 
DBH (see Table 1) for the plot ‘Pine 1’ - 24,4 cm (max 
65, min 5), for ‘Pine 2’ - 25,9 cm (max 60, min 7) and 
28,3 cm (max 56, min 6 cm) for the plot ‘Pine 3’. 

All the research areas were facing eastwards and 
the ground slope was 15-17°. All the species observed, 
originated from self-sowing, except beech located in the 
‘Pine 1 ’ plot. Beech in that area originated partially from 
sowing by human. 

Statistical analyses were performed using 
STATISTICA software. In all three areas, DBF! values 
are coming to form a normal, right-skewed distribution 
(Shapiro-Wilk test). 

Each research area was a 25 x 25 m rectangle. Such 
area (625 m 2 ) is appropriate for gap models (Bugmann 
1994). Running 200 simulations for each research area 
is equivalent to analyzing an area of approximately 12,5 
ha (200x0,0625). 

Each research area was constructed using standard 
measuring tape to establish a rectangular coordinate 
system (X axis along E-W and Y axis along N-S line). 
Each and every tree within such established research 
area was assigned a unique number and measure 
(DBH, height, crown radius, age). The position of 
every tree was also established at 5 cm accuracy and 
then marked on a map using orthogonal survey method. 
The age was established with the help of Pressler drill. 
The diameter at breast height was measured using 
a dedicated meter at 1 ,3 m over the ground level. The tree 
height was measured with Leiss BL8 height-meter. In 
order to verify leaf area index (LAI) calculations, three 
hemispherical photographs were taken per research 
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Table 1 . Characteristic of trees on research plot 



Plot 


Species 


D [cm] 


H [m] 


Age (years) 


N 


mean 


min 


max 


mean 


min 


max 


mean 


min 


max 


1 


Abies alba 


10,0 


5,0 


18,0 


7,6 


4,0 


13,0 


14,8 


11,0 


16,0 


8 


1 


Fagus sylvatica 


8,3 


6,0 


14,0 


5,4 


4,0 


9,0 


15,4 


12,0 


20,0 


7 


1 


Betula pendula 


20 


20 


20 


16 


16 


16 


26 


26 


26 


1 


1 


Picea abies 


9,4 


4,0 


21,0 


8,6 


4,0 


18,0 


16,2 


11,0 


35,0 


10 


1 


Pinus sylvestris 


24,4 


5,0 


65,0 


18,0 


5,0 


39,0 


67,5 


10,0 


148,0 


35 


2 


Abies alba 


7,0 


7,0 


7,0 


7,0 


7,0 


77,0 


19,0 


19,0 


19,0 


1 


2 


Fagus sylvatica 


7,0 


5,0 


8,0 


5,2 


4,0 


6,0 


14,0 


4,0 


6,0 


5 


2 


Picea abies 


12,4 


9,0 


16,0 


10,0 


8,0 


12,0 


33,2 


24,0 


43,0 


5 


2 


Pinus sylvestris 


25,9 


7,0 


60,0 


14,9 


5,0 


29,0 


68,9 


19,0 


160,0 


37 


3 


Abies alba 


5,0 


5,0 


5,0 


6,0 


6,0 


6,0 


19,0 


19,0 


19,0 


1 


3 


Fagus sylvatica 


7,0 


5,0 


8,0 


5,2 


4,0 


6,0 


14,0 


10,0 


16,0 


5 


3 


Picea abies 


7,6 


5,0 


13,0 


6,8 


5,0 


10,0 


20,2 


13,0 


35,0 


5 


3 


Pinus sylvestris 


28,3 


6,0 


56,0 


16,1 


7,0 


28,0 


75,4 


19,0 


149,0 


39 



Designations: D - d.b.h., H- height, N - number of trees 



area (model FORKOME, SLR camera Canon EOS 
5D 12MP, Sigma 8mm f/3,5 DG FISH EYE lens with 
180° viewing angle). The hemispherical photographs 
were analysed with the help of the Gap Light Analyzer 
software (Frazer et al. 2000). The software enabled 
measuring not only LAI, but also the absolute amount of 
direct and scattered light, structure of canopy openness, 
spatial variability of canopy foliage and biomass of the 
assimilation apparatus. Next, the obtained results were 
used to verify the FORKOME model input parameters 
(assimilation apparatus biomass). Selected data were 
stored in csv format (using Microsoft Excel spreadsheet) 
and supplied to the FORKOME model. 

The FORKOME model is based upon the theory of 
forest ecosystems patch dynamics. It can be used to simulate 
forest association succession while taking into account 
the association characteristics (species composition, age, 
etc.), as well as the specificity of individual growth. The 
latter can be simulated as a function of both external 
(light intensity, ambient temperature, etc.) and internal 
(competition) ecological factors. 

Based on appropriate equations, the model simulates 
reproduction, growth and death of trees every year 



(Brzeziecki 1999; Kozak et al. 2012). For every species 
available in the model, specific death age has been 
assumed, e.g. 400 years for pine. The sets of parameters 
used in the model take into account the most important 
ecological and breeding properties of each species. 
Since the parameters reflect several properties and 
habitats of specie, they can vary to a certain degree. The 
FORKOME model is also capable of taking into account 
natural disasters, forest fires in particular (Kozak et al. 
2012). In general, pine requires light for successful 
restoration (weak natural renewal under forest canopy 
and optimal on fire -opened spaces). The model allows 
simulating both types of restoration. It is possible to 
display the distribution of trees in the research areas 
as shown in Figure la, lb, lc. Simulation time was 
set to 100 years and simulation type to Monte Carlo 
(200 runs). Forecasts of biomass changes and the 
number of trees were made for all pine areas. 

To analyze climate changes, the FORKOME model 
uses the method proposed by Kahn (1994) which has 
been used previously, e.g. in the SILVA hybrid model of 
forest dynamics (Pretzsch et al. 2002). Out of the factors 
considered by Kahn (1994), only the basic ones were 
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Figure 1. Beginning stage (3D visualization in FORKOME model) of plots in control scenario: a - ‘Pine 1 ’, b - ‘Pine 2’, c - ‘Pine 3’ 



used in this study: annual temperature amplitude, mean 
vegetation period temperature, and vegetation period 
total precipitation. 

The influence of climate changes (total effective 
temperatures and total annual precipitation) on the rate 
of pine replacement by other species (beech, fir, spruce, 
hornbeam) was also the subject of this study. To this 
end, ‘temperature changes scenario’ and ‘precipitation 
changes scenario’ tools were used. To assess the influence 
of temperature changes the following values were used: 

the plot ‘Pine 1’: control - 1425 degree days, ‘warm- 
dry’ and ‘warm-humid’ scenario - 1625 degree days, 
‘cold-dry’ and ‘cold-humid’ scenario - 1225 degree 
days. 

- The plot ‘Pine 2’: control - 1435 degree days, 
‘warm-dry’ and ‘warm-humid’ scenario - 1635 degree 
days, ‘cold-dry’ and ‘cold-humid’ scenario - 1235 
degree days. 

- The plot ‘Pine 3’: control - 1440 degree days, 
‘warm-dry’ and ‘warm-humid’ scenario - 1640 degree 
days, ‘cold-dry’ and ‘cold-humid’ scenario - 1240 
degree days. 

Similarly, to assess the influence of precipitation 
changes, the following values were used: 

- The plot ‘Pine 1’: control - 718 mm (data from 
weather station in Zwierzyniec, Poland), ‘warm-dry’ 
and ‘cold-dry’ scenario - 618 mm, ‘warm-humid’ and 
‘cold-humid’ scenario - 818 mm. 

- The plot ‘Pine 2’: control - 713 mm (data from 
weather station in Rava-Rus’ka, Ukraine), ‘warm-dry’ 
and ‘cold-dry’ scenario - 613 mm, ‘warm-humid’ and 
‘cold-humid’ scenario - 813 mm. 



- The plot ‘Pine 3’: control - 645 mm (data from 
weather station in Lviv, Ukraine), ‘warm-dry’ and 
‘cold-dry’ scenario - 545 mm, ‘warm-humid’ and ‘cold- 
humid’ scenario - 745 mm. 

No climate change scenario (kn) was used as 
control. Next, the following climate change scenarios 
are defined: 

- warm-dry (wd): temperature increase by 200 
degree days and decrease in precipitation by 100 mm, 

- warm-humid (wh): temperature increase by 200 
degree days and increase in precipitation by 100 mm, 

- cold-dry (cd): temperature decrease by 200 
degree days and decrease in precipitation by 100 mm, 

- cold-humid (ch): temperature decrease by 200 
degree days and increase in precipitation by 100 mm, 

3. Results 

Diagrams presenting results of prognosis of plots 
‘Pine 1’, ‘Pine 2’, and ‘Pine 3’ are placed in the paper 
under the control scenario, as well as in 4 climate 
changes scenarios. Analyses of the aforementioned 
diagrams indicate that under the control scenario, (cs) 
pine biomass will increase during the period predicated 
at 56 (Fig. 2a) or 60 years (Fig. 3a, Fig. 4a). In ‘Pine 
1’, the area placed west, pine biomass will be 420 t/ha, 
and in the other plots, i.e. ‘Pine 2’ and ‘Pine 3’, pine 
biomass will increase to 550 and 670 t/ha. Decreasing 
tendency of pine biomass may be observed in almost all 
areas under the control scenario after the 60th year of 
the prognosis. Flowever, inconsiderable increase of pine 
biomass occurs at the end of prognosis time (between 






I. Kozak et al. / Forest Research Papers, 2013, Vol. 74 (3): 215-226. 



219 



Figure 2. Biomass total production, 
accumulation (a) and number 
of trees (b) of analysed species 
in the control scenario on the plot 
‘Pine 1’ 




Figure 3. Biomass total production, 
accumulation (a) and number 
of trees (b) of analysed species 
in the control scenario on the plot 
‘Pine 2’ 




Figure 4. Biomass total production, 
accumulation (a) and number 
of trees (b) of analysed species 
in the control scenario on the plot 
‘Pine 3’ 
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Figure 5. Biomass total production, 
accumulation (left) and number 
of trees (right) of analysed species 
in the warm dry climate scenario 
(wd) on the plots (a, b) ‘Pine 1 
(c, d) ‘Pine 2’, (e, f) ‘Pine 3’ 
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Figure 6. Biomass total production, 
accumulation (left) and number 
of trees (right) of analysed species 
in the warn humid climate scenario 
(wh) on the plots (a, b) ‘Pine 1’, 

(c, d) ‘Pine 2’, (e, f) ‘Pine 3’ 
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Figure 7. Biomass total production, 
accumulation (left) and number 
of trees (right) of analysed species 
in the cold dry climate scenario (cd) 
on the plots (a, b) ‘Pine 1’, (c, d) 
‘Pine 2’, (e, f) ‘Pine 3’ 
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Figure 8. Biomass total production, 
accumulation (left) and number 
of trees (right) of analysed species 
in the cold humid climate scenario 
(ch) on the plots (a, h) ‘Pine 1’, 

(c, d) ‘Pine 2’, (e, f) ‘Pine 3’ 
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scenarios. Changing climate scenario, the cold dry 
climate scenario (cd) indicates that pine biomass will 
come to 450 t/ha in the plot ‘Pine 1’, and 690 t/ha and 
670 t/ha in the plots ‘Pine 2’ and ‘Pine 3’. Beech and fir 
biomass will decrease, but spruce biomass will increase. 
Regarding the number of trees, spruce will increase in 
comparison to the control, especially, under the plot 
‘Pine 1’ (Fig. 7b) and under the plot ‘Pine 2’ (Fig. 7d). 
Number of beech trees will remain at a minimum level. 

In the cold humid (eh) climate scenario, pine 
biomass will increase in comparison to the control. It 
will be even larger compared to the cold dry climate 
(cd) scenario. It will come to 520 t/ha in the plot ‘Pine 
1’ (Fig. 8a), 700 t/ha in the plot ‘Pine 2’ (Fig. 8c), and 
670 t/ha in the plot ‘Pine 3’ (Fig. 8e). Regarding number 
of trees, spruce increase can be observed. Beech and fir 
individuals will remain at a low level. 

4. Discussion 

We carried out simulations with the FORKOME 
model in 3 research areas of 25x25 m (1/16 ha) each. 
Therefore, within 625 m 2 of the area, 200 used simulation 
sequences correspond with equilibrium landscape of 
12,5 ha. 

In other types of models, for instance the JABOWA 
model (Botkin et al. 1972), 10x 10 m areas were used with 
100 simulations. The smallest area was 0,5 ha and was 
used in the FORSKA model (Leemans, Prentice 1987). 
Research carried out in Roztocze proved that it is relevant 
to use research area 625 m 2 of size in patchwise models. 

Results of studies carried out with the FORKOME 
model indicate that pine will not be dominated by 
other species during the period of the next 100 years. 
Studies presented in literature have corroborated 
the aforementioned thesis. For instance, it has been 
observed that pine in the area close to the South border 
of pine specie range, in Northeast Spain (Martinez- 
Vilalta et al. 2008), increased the sectional area by 84% 
during the 20th century. It has not been excluded that 
the trend relates to increasing of CO, in the atmosphere. 

Simulation of the dynamics of pine stands in 
Roztocze was based on research areas placed in pine 
stands where podzol soil dominates. This type of soil 
is usually overgrown by pine stands (for instance in 
Roztocze), while spontaneous combustion of the forest 
enable pine to restore (Kozak et al. 2012). 

It is proved from computer simulations that the 
number of pine trees decreases continually - old trees 



whither, and canopy cover does not allow younger 
individuals to receive enough amount of light for 
growing. Diagrams indicate slow tree restoration 
(without a fire). 

Other authors have often used the FORKOME model 
in computer simulations regarding different climate 
change scenarios (warm dry, warm humid, cold dry and 
cold humid) (Brzeziecki 1999). It has not only been 
a theoretical application, but also practical application, 
because it enables determining the directions of changes 
in the environment and then, defines proper management 
and protection strategies. 

It is not easy to predict climate changes for a long 
period of time, because it is impossible to assume that 
climate will change only in one direction. That’s why 
we carried out four simulation variants and results 
are convergent with the pine ecological requirements 
(Ellenberg 1986; Kahn 1994). 

Pine is characterised by broadest ecological 
spectrum among all analysed species. It is important in 
particular in case of extreme temperatures and humid 
conditions in connection with weak podzol soil. Beech, 
fir, and spruce die as first in such conditions. Prognosis 
results done with the FORKOME model are relatively 
representative. Studies carried out in Roztocze 
National Park (Maciejewski 2009; Maciejewski 
2011; Maciejewski, Szwagrzyk 2011), Rava-Rus’ka 
Landscape Reserve and (Cherniavskyi, Sawka 2004) 
and Yavorivskyi National Park (Horoszko, Flomiuk 
2011) can corroborate these results. 

Studies with the FORKOME model give results 
which are convergent with results presented in literature 
(Maciejewski 2009; Maciejewski 2011; Maciejewski, 
Szwagrzyk 2011; Horoszko, Homiuk 2011). During the 
period of 100 years, the FORKOME model predicts that 
pine stands will remain, even in the warm dry climate 
scenario (wd), and what is corroborated by the pine 
classification, as the specie is relatively adapted for 
drought (Ellenberg 1986). It predicts that changes in 
species’ structure of pine stands will occur by increasing 
number of beech, fit, and hornbeam in warming climate 
scenarios. It is worth mentioning that warming climate 
process is one of hypothesis (Fabijanowski, Jaworski 
1996) to explain that phenomenon. We made an 
assumption that climate conditions for analysed pine 
stand areas in Roztocze in the FORKOME model are 
compatible with conditions in the research areas. These 
assumptions proved to be propitious for pine. 

Prediction research on stands development, in 
regard to climate change conditions, is a significant type 
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of scientific research. The FORKOME model is a good 
enabler to obtain results that prove the thesis assumed in 
our paper. Results obtained from the prediction research 
correspond to results obtained from field works and 
phenomenon described in literature. The FORKOME 
model proved to be efficient tools adjusted to the 
research regarding dependences and changes in species 
structure of forest stands, species biomass, and number 
of trees. It is possible to devise methods in the model 
to analyze many parameters of forest stands conditions 
at the same time what allows obtaining precise results. 
The FORKOME model enables analyzing the dynamics 
of forest stands in regard to climate change conditions 
during any period of time, chosen by the user of the 
model. 



5. Conclusions 

Computer simulations presented in the paper indicate 
that pine will be basic species regardless of climate 
changes, and in case of worse soil and water conditions, 
it will be the main specie of forest stands grown on 
podzol soil in Polish and Ukrainian parts of Roztocze. 
Pine will not become less competitive in Roztocze in 
comparison to beech or fir, especially in stands grown 
on weak podzol soil. Beech, fir, and spruce grown in 
these stands will be always second-rated and considered 
admixture species. 

In pine stands in Roztocze, pine will maintain its 
domination regarding biomass in the control scenario, 
warm dry climate scenario (wd), warm humid climate 
scenario (wh), cold dry climate scenario (cd), and cold 
humid climate scenario (eh) in the period of the next 
100 years. Pine has optimal growing conditions in 
Polish and Ukrainian parts of Roztocze. 

Thanks to adequate tools (climate change scenarios), 
the FORKOME model is adjusted to prediction research 
regarding forest studies and it enables verifying assumed 
hypothesis. Results obtained from prediction research 
made in the FORKOME model proved the usefulness of 
the model for simulation of dynamics of pine stands in 
Polish and Ukrainian parts of Roztocze. 
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